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Abstract 

We use the leading singularity technique to determine the planar six-particle two-loop MHV 
amplitude in 7\A = 4 super Yang-Mills in terms of a simple basis of inte gra ls. Our result for the 
parity even part of the amplitude agrees with the one recently presented in [18|. The parity-odd part 
of the amplitude is a new result. The leading singularity technique reduces the determination of the 
amplitude to finding the solution to a system of linear equations. The system of equations is easily 
found by computing residues. We present the complete system of equations which determines the 
whole amplitude, and solve the two-by-two blocks analytically. Larger blocks are solved numerically 
in order to test the ABDK/BDS iterative structure. 
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I. INTRODUCTION 



Scattering amplitudes in massless gauge theories are remarkable objects with many prop- 
erties hidden in the complexity of their Feynman diagram expansion. It has been known for 
decades that much information about an amplitude can be gleaned just from knowing the 
structure of its singularities (see jl|). In Yang-Mills theories, extensive use of branch cut 



singularities was shown to tame much complexity in the computation of loop arnplitudes via 
techniques that came to be known as the unitarity based method [2], [s, Q, 0, ^ Q, [s]. 

One of the most surprising features of Yang-Mills theory and gravity is that their tree 
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y determined by exploiting only their behavior near a 



111 ]. Another surprise, this time in A/" = 4 SYM, is that 



level amplitudes can be cornp 
subset of their singularities [q, 
the problem of computing one-loop amplitudes can be reduced to that of computing tree 
amplitudes [l2|. The key to such striking simplification is that although the loop amplitude 
possesses poles and many branch cuts with a complicated structure of intersections, it is 
completely determined by the structure of the highest codimension singularities. These are 
known as the leading singularities 1| and are computed by cutting all propagators in the 
diagram. 

Applying the same technique at higher loops was first attempted in [l^ and refined 
15|. In 13] and in |l^ both the leading singularity as well as subleading {i.e., lower 
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codimension) singularities were used to constrain the form of higher loop amplitudes. In 16| , 
the leading singularity was shown to be much more powerful than expected. It turns out 
that in massless theories, whenever all propagators are cut, one is actually studying many 
iso 



m 



ated singularities at the same time. The proposal of [16|], building on observations made 
isl . is to use each isolated singularity independently 



has three notable features. Firstly, 



The new leading singularity technique outlined in 
for any amplitude under consideration it builds a natural set of integrals, which we call 
the geometric integrals, that can be used to construct a basis. Secondly, the coefficients of 
the integrals can be determined by solving linear equations. Finally, these linear equations 
a.e easily obtained by computin^residues using Cauchy's theorem. The utility of this uew 
technique was demonstrated in [16'], where it was shown that it easily reproduces the result 
for the two-loop five-particle amplitude in A/" = 4 SYM previously computed in 
the unitarity based method. 
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In this paper we apply the leading singularity technique to a much more challenging 

(2) 

case, the planar two-loop six-particle MHV amplitude AqI^-^y. The parity-even part of 
the normalized amplitude ^g^MHvMe^MHV computed recently in 1&\ using the unitarity 
based method. This was already an impressive display of computational power. In the 



five particle case studied in 
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24| the parity-odd part is noticeably of a higher degree of 



complexity than the parity-even part, and there is no reason to suspect that this would not 
be the case also for n = 6. Here we find that the full coefficients (both the even and odd 
parts) emerge by solving the relatively simple linear equations presented explicitly below. 
For six particles we find a new phenomenon which was not encountered in the n = 4, 5 

n 

cases studied in [16|: while there are 177 geometric integrals, the equations only fix 159 
linear combinations of their coefficients, leaving 18 linear combinations undetermined. The 
reason is that the set of geometric integrals is overcomplete, so we cannot have expected 



to find a unique solution. One can show that by using well-known techniques [19|, |20| it is 
possible to build linear relations between seemingly independent integrals. In section IV we 
analyze the relevant reduction identities and identify 18 relations amongst the integrals in 
the geometric basis. Thus there is no ambiguity beyond that required by reduction identities. 



so we conclude that A, 



(2) 
,MHV 



is in fact completely determined by its leading singularities. 



It is important to stress that the leading singularity method turns loop integrals into 
contour integrals which are finite in four dimensions and knows nothing about how one might 
choose to regulate the infrared divergences that typically appear when carrying out the loop 
integrals. In dimensional regularization, amplitudes occasionally contain additional terms 
in the integrand which vanish in Z) = 4. These so-called "//-terms" (see for example j21| for 
a thorough treatment) cannot be detected by the leading singularity in D = 4. 

One motivation for computing the MHV six-particle amplitude, beyond its serving as a 
testing ground for the leading singularity method, is to study the proposed iterative relation 
between MHV loop amplitudes known as the ABDK/BDS ansatz 22, l23|]. The ansatz has 
been shown to hold for four particles up to three loops 22, l23|] and for five particles up to 



two loops 171, l2J]. However, it was shown to break down for the parity-even part of the 



two-loop six-particle MHV amplitude in 



18|. In this paper we find numerical evidence that 



the parity-odd part of the amplitude does satisfy the ABDK/BDS ansatz. 
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II. OUTLINE OF THE CALCULATION 



The object of interest is the planar six-particle two-loop MHV amplitude in = 4 
super Yang-Mills. The goal is to find a compact expression for this amplitude as a linear 
combination of relatively simple integrals. The leading singularity method [l^ provides 
both a natural set of integrals to work with, as well as a system of linear equations which 
determine the coefficients of those integrals. 

In this section we provide a detailed outline of the steps involved in setting up the calcu- 
lation. The first three subsections are relevant to NMHV as well as MHV amplitudes, since 
the homogeneous part of the system of linear equations is helicity independent. In subsec- 
tion II. D we compute the inhomogeneous terms for the MHV helicity configuration. The 
final linear equations which determine the coefficients of the MHV amplitude are presented 
explicitly in section HI. 

A. Review of the Leading Singularity Method 

Suppose we are interested in calculating some L-loop scattering amplitude A. On the 
one hand, the amplitude may of course be represented as a sum over Feynman diagrams Fj, 



where k are external momenta and ia are the loop momenta. However it is frequently the 
case, especially in theories as rich as A/" = 4 SYM, that directly calculating the sum over 
Feynman diagrams would be impractical. Rather the calculation proceeds by expressing A 
as a linear combination of relatively simple integrals in some appropriate basis {/«}, 



and then determining the coefficients Cj by other means, such as the unitarity based method. 
If the set of integrals {/j} is overcomplete, then the coefficients Ci{k) are not uniquely defined. 

The basic idea underlying the leading singularity method is that the sum over Feyn- 
man diagrams in (12. ip possesses singularities which must be properly reproduced by any 
representation (12. 2p of the amplitude in terms of simpler integrals. At the same time, any 




(2.1) 




(2.2) 
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singularities in the set of integrals which are not present in the sum over Feynman diagrams 
must be spurious. 

The most common kind of singularities in Feynman diagrams are poles, associated to 
coUinear or multi-particle singularities, and branch cuts, associated to unitarity cuts. These 
branch cuts can themselves possess branch cuts leading to higher codimension singularities. 



The latter are computed by cutting propagators or equivalently [l5| by promoting the loop 



integral to be a contour integral. (The observation that the Lorentz invariant phase space 



integral of a null vector can be recast as a contour integral was first discussed in [25|.) The 
contour is chosen to reproduce the behavior of the delta-functions in the cut calculations. In 
general, this gives rise to contours which compute the residue on several isolated singularities 
at the same time. 

For example, consider the one-loop massless scalar box. Replacing each propagator by a 
delta-function leads us to consider the integral 



d^i 6i£')5i{£ - hy)6ii£ -h- hy)6iii + hY) . (2.3) 

For generic external momenta ki these delta-functions localize the i integral onto two discrete 
points in complex £-space (C^). With the leading singularity method we do not replace 
propagators by delta-functions, but rather we consider two separate contours in C^, each 
of which computes the residue of the integrand on only one of the two isolated singularities. 
Then by equating (12. ip and (12.21) and performing the integral 

Vq(A;) [ dHh{k,e)= [ dHy^F.ikJ) (2.4) 

we obtain one linear equation on the coefficients Cj for each contour F. At L loops each 
contour is a T^^ inside C^^. Since the number of isolated singularities in a generic L- 
loop diagram can be as high as 2^ (simple diagrams can have fewer isolated singularities), 
the leading singularity method gives rise to an exponentially large (in L) number of linear 
equations for the coefficients q. We note that the homogeneous part of these linear equations 
(the left-hand side of f l2.4p ) depends only on the set of integrals {/j} and the choice of 
contours, while the details of which particular helicity configuration is under consideration 
enters only into the inhomogeneous terms on the right-hand side. 
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FIG. 1: The five indepedent 8-propagator topologies. Each diagram represents a sum of those 
Feynman diagrams in which all of the 8 indicated propagators are present. In each diagram the 
external momenta are labeled clockwise beginning with ki at the position of the arrow. Also in 
each diagram p is the loop momentum in the left loop and q is the momentum in the right loop. 

B. Choosing Useful Contours 

The formula (12. 4p gives a linear equation on the coefficients c, for any contour T in C^^. 
Of course if we choose some random contour T then we will typically get the useless equation 
= 0. In order to get useful equations we should use T^-^ contours which calculate residues 
at the known singularities of the right-hand side. It is clear that in the sum over Feynman 
diagrams, singularities occur when internal propagators go on-shell. 

For six particles at two-loops there are two classes of useful contours. The most 
obvious contours are those which are chosen to calculate the residue at points in 
where eight propagators go on-shell simultaneously. These contours are associated with the 
five different topologies shown in fig. [H Actually each topology in fig. [1] is a diagrammatic 
shorthand for four distinct contours. For example, the singularities of Feynman diagrams 
with topology (D) are situated at the locus 

S(D) = {{p, g) G X e : p2 = 0, {p + hf = 0, {p + = 0, (p - k^f = 0, 

g2 = 0, {q + kif = 0, {q + ki2f = 0, {p + qf = 0} . (2.5) 

For generic external momenta ki S(^o) consists of four distinct points in of the form 
(p(*)^ g(i)) for i^j = 1,2. Correspondingly there are four different contours F associated with 
topology (D), one which computes the residue of the integrand at each of these four isolated 
singularities. 

The less obvious contours are those in which only seven propagators are apparent but 
an eighth singularity appears due to a Jacobian. These contours are associated with the 
eight different topologies shown in fig. O For example, let us consider topology {F). For 
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FIG. 2: The eight independent 7-propagator topologies. See fig. [T] for details. 

fixed loop momentum p the singularities in the q integral occur at the locus 

S^F), = {g G : = 0, (g + = 0, (g + ^2? = 0, (g + p)^ = 0} , (2.6) 

which consists of two points {g^^^g^^^} in C'^. For each of these two singularities there 
is a contour Vq such that integrating g over Vq computes the residue at that singularity. 
Integrating over either contour produces the same Jacobian factor 

1 1 1 



(2.7) 



g2(g + k^f^q + k^^f^q + pf 2 {k^ + k^f^p - k^f " 
The new singularity l/{p — kiY combines with the three remaining singularities manifest in 
topology {F) so that the integral over p can be localized by integrating over contours which 
compute the residue at the points 

= G : p2 = 0, (p + k^f = 0, (p - k^2f = 0, (p - k^f = 0} . (2.8) 

We proceed analagously for each of the eight topologies shown in fig. [21 In each case 
we first integrate the right-hand loop momentum g and then use the additional singularity 
generated by the Jacobian to integrate the left-hand loop momentum p, thereby completely 
localizing the integral onto a set of discrete points. 

For the MHV amplitude it turns out that the linear equations generated by the 13 types 
of contours described in figs. [U and [2] are sufficient to uniquely determine all coefficients, 
although there certainly are additional contours that could be used to generate additional 
equations from (12.41) . We have checked that a class of additional equations are indeed sat- 
isfied by the MHV coefficients, thus providing a strong consistency check on the coefficients 
obtained from solving the equations in section III. 
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FIG. 3: The geometric integrals for the two-foop six-particle amplitude in A/" = 4 SYM. In each 
diagram the external momenta are labeled clockwise beginning with ki at the position of the arrow 
and the number in parentheses denotes the number of independent permutations of the diagrams. 
As discussed in section IV, this is an overcomplete set: the 177 integrals here obey 18 linear 
relations. 

C. The Geometric Integrals 



Although equation fl2.4l) can be used to generate linear equations for coefficients in an 
arbitrary basis {/i}, the leading singularity method suggests a natural set of integrals in 
which the left-hand side is easy to compute and the individual integrals have a geometric 
interpretation, in a sense we now explain. 

The procedure to determine the natural set of integrals in which to represent an amplitude 
starts by realizing that in A/" = 4 SYM tadpoles, bubbles, and triangles are unneccesary 
(see I^gI] for a thorough discussion). This means that we have to start by considering all 
topologies of sums over Feynman diagram with no triangles or bubbles. At one-loop, only 
sums of Feynman diagrams with the topology of a box are needed. For six particles at two- 
loops, we find five topologies with eight propagators, shown in fig. [H and eight topologies 
with seven propagators, shown in fig. [21 

After each topology is identified, a first approximation for reproducing all of the leading 
singularities is to use just the scalar integrals with all of the appropriate topologies. In 
general, it turns out that such integrals are not enough to reproduce the singularities of 
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Feynman diagrams. This will manifest itself in the failure of the linear equations (12 ■4p to 
have any solutions, indicating that the set of integrals must be enlarged. 

The next step is to introduce scalar integrals with additional propagators. At one-loop 
this step gives rise to pentagons in addition to boxes, which turns out to be sufficient for 
any n. At two loops we add the scalar pentagon-pentagon integrals shown in fig. [31 At this 
stage some of the equations are solved {i.e. some of the leading singularities are correctly 
reproduced), while others are not. 

It is then necessary to supplement additional integrals which must have non-zero residue 
on the missing singularities and zero residue on the ones which already work, in order to 
avoid spoiling them. The way to ensure that one has zero residue on a given pole is to 
include a zero in the numerator of the integrand which cancels the corresponding pole. In 
this form integrals with scalar numerators appear. Note that only numerators which cancel 
poles appear naturally. 

The process of expanding the set of integrals by including additional numerator factors 
ends when one is able to solve all of the equations (12. 4p . We call the integrals that are 
naturally constructed in this manner geometric integrals. The set of geometric integrals 
might be overcomplete if the equations do not determine a unique solution. For the six- 
particle MHV amplitude at two loops this process leads to the 177 geometric integrals shown 
in fig. [3l For NMHV amplitudes it is possible that additional integrals, such as a pentagon- 
boxes with two numerators, might be required. 

D. The Inhomogeneous Terms for the MHV Amplitude 

As mentioned above, it is obvious from (12. 4p that the homogeneous part of the system of 
linear equations is helicity independent, so the above discussion of the contours and geometric 
integrals applies to both the MHV and NMHV configurations. The inhomogeneous part on 
the right-hand side of (12. 4p is easily obtained for any contour F by computing the product 
of tree amplitudes sitting at the 'blobs' in the corresponding topology in fig. [1] or fig. [2l 
This product is evaluated at the value {p^^\q^^^) of the loop momenta where the contour F 
localizes the integral. 

For MHV configurations there is an enormous simplification since it turns out that this 
product of tree amplitudes always comes out to be or 1 times the corresponding tree-level 
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amplitude, as can easily be shown by using the technique introduced in ISj where sums over 



the full J\f = 4 supermultiplet in the internal lines, which complicate the computation 
Q, are automatically done by using simple identities. 

For the contours associated with the 13 topologies shown in figs. [T] and [2] we find that the 
right-hand side of fl2.4p for an MHV configuration are 

(A) : 5(p,i)S{qA)A6MBV^ 

(B) : '^(p,6)5(g,l>^6'MHV' 

(C) : 0, 

{D) : '^(p,6)^6^MHV ' 

{E) : '^(p,6)'^(g,2>^6'MHV) 

(F) : ^{p,6)^6^}iY ' 

(G) : 0, 

(H) : 0, 

{!) '■ '^(p,6>^(p,g>^6^MHV ' 

(J): 0, 

(K) : 5[p,6Aq,i]AQ^Mnv , 

(L): 0, 

(M) : 0. (2.9) 

The Kronecker deltas such as 5(p,i) for topology (A) arise because, in that example, the 
solution for p is either of the form Ap oc Ai or Ap oc Ai, and the product of tree amplitudes 
vanishes on latter solution. 



III. THE MHV EQUATIONS AND THEIR SOLUTIONS 

We now assemble all of the ingredients prepared in section II for the MHV amplitude. 
By evaluating ( 12. 4p on all of the contours associated with the topologies shown in figs. [T] 
and [2] (together with all of their cyclic and mirror- image permutations), and using (12.91) on 
the right-hand side, we find a system of 396 linear equations for the 177 coefficients of the 
geometric integrals in fig. [31 

Generically, 396 linear equations in 177 variables have no solution, but in this case we find 
that the equations are in fact underdetermined: they only fix 159 linear combinations of the 
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177 coefficients, leaving 18 free parameters. In other words, we find that there are 18 hnear 
combinations of the geometric integrals in fig. [3] which have vanishing leading singularity 
on all of the contours described by figs. [T] and [2J One logically possible conclusion from 
such a result might have been that the leading singularity method is not enough to uniquely 
determine the two-loop six-particle MHV amplitude, which would have been disappointing. 

Fortunately, as mentioned in the introduction, it turns out that integral reduction iden- 
tities imply that the set of geometric integrals is overcomplete. In other words, there are 
linear combinations of the 177 geometric integrals which not only have vanishing leading 
singularity but actually vanish identically. In section IV we analyze these reduction identi- 
ties and show that there are 18 linear combinations of the geometric integrals which vanish, 
precisely accounting for the abovementioned ambiguity in solving the leading singularity 
equations. The conclusion is therefore that the two-loop six-particle MHV amplitude is in 
fact uniquely determined by knowledge of its leading singularities. 



A. Presentation of the Equations 

In order to demonstrate that the leading singularity method is not just black magic, we 
present here explicit expressions for the equations which determine all coefficients mhv' 
including both the parity-even and odd terms. As just discussed, 18 coefficients out of 177 
are actually redundant. In order to somewhat simplify the presentation of the equations we 
will choose a convenient 'gauge' which uniquely fixes all of the ambiguity. This amounts to 
choosing a basis of geometric integrals. 

Several such choices are possible; the choice we make here is to spend the 18 gauge 
parameters by setting to zero the 6 coefficients Ci2 and the 12 coefficients C24. Once this is 
done all remaining coefficients are uniquely determined. A nice advantage of this choice is 
that several other coefficients turn out to also vanish identically: the 12 coefficients C22, the 
12 coefficients cis, the 6 coefficients cie and the 3 coefficients C23 are all zero. 

Ultimately then there are only 126 nonzero coefficients, associated with just 15 out of 
the 21 integrals shown in fig. [3l We label the j-th permutation of coefficient Cj as c^p . 
Permutation j maps the labeling (1, 2, 3, 4, 5, 6) of the external momenta to: 

1 : (1, 2, 3, 4,5,6), 2 : (2, 3, 4, 5,6,1), 3 : (3, 4, 5, 6,1,2), 4 : (4, 5, 6, 1, 2, 3) , 
5 : (5, 6, 1, 2,3,4), 6 : (6, 1,2,3,4,5), 7 : (6, 5, 4, 3,2,1), 8 : (5, 4, 3,2,1,6), 
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9 : (4,3,2,1,6,5) , 10 : (3, 2, 1, 6, 5, 4) , 11 : (2, 1, 6, 5, 4, 3) , 12 : (1, 6, 5, 4, 3, 2) (3.1) 



Since all of the coefficients for the MHV amplitude are proportional to the tree amplitude, 
we can go ahead and divide the right-hand side of all equations by ^q'^mhv Equivalently we 
can say that solving the equations below yields the integral coefficients for the normalized 
amplitude ^g^MHvMe'^MHV ^ ^^^^ comment is that we move all of the Jacobian factors (see 
for example eq. (12.71) ) to the right-hand side of the equations. 

Finally we are ready to present the equations obtained by considering the contours asso- 
ciated with the topologies in figs. [1] and [2l In each equation p and q are understood to be 
evaluated at the locations {p^^\q^^^) of all the leading singularities. 
Topology A: 

(2) (P - fc234)'(g - fc34)'cg + (p + kei)^{q - hufcf^ _ ^ ^ . . 

Ci + {p + q- k2Zif ~ 4Si2S23S45S56')(p,l>')(g,4> (3.2) 

Cl9^ + (P - kifcf = Asuip - kiY{si23S3A5 " Si2S45)S{pfi) (3.3) 



Topology D: 



Topology E: 



(3) , . , x2 (3) , jP - ki2)\q + fc234)^cS3^ + (p " ki)\q - fceO^cg 
+[p- fci2j + (q-k^^ 



Topology F: 



Topology H: 



1) 

4S23S45S56(P - ^12)''5(p,6>5(g,2> (3.4) 



cf + ^^^^4.W^(.a) (3.5) 



(6) , (p-fcl2)^(9-fc56)M? + (p-fcl)^(g-fc6)Mg 

{p + k,,f{q + k,2f 
cf^^{q-k,fcf , 4^^ + (p-fc0^cf) _ 

^ {q + kuy + {p+k,ey ~ 



Topology I: 



.(6) 



{p - ku)\q - k,,fcf^ + {p- k,)\q - k,Ycf^ 



^ {p + k,enq + k,y 

, 4g + iq- ke?cf^ 4g + {p- kr^ <^^ 
{q + k^f ^ {p + kr^^^ 



f^^h )2 = 4si23(si23S345 " Si2S45)(5(p,6><5(p,g> (3.7) 

^ ^56 j 
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Topology J: 



, , , , X2J12) (5) / , x2 (5) 

(5) + l^p - fceij c^o C21 + (^g - ftsej c^q 

C7 n 7 TTT^ r 
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(p-A;6)2(g + fci23)2 



+ 



{p - keifiq - k45efc^r] + {p - ke)^{q - k^efc^^ ^ ^ 



Topology K: 



{p + ku,nq-ker 



^(6) ^ +{P - kufjq - A:56)^4f + {p - k,f{q - k.fcfj 
{p + ki56y{q + ku3y 

■ + ^° + t , 1 - 4.,2.56(? - ^6)^5(,,,) (3.9) 



{q + ki23Y (P + ^456)^ 



Topology L: 



(3) , (p-M^(g-fc56)M? + b-fel)^(g-fc6)Mg 



(p + /i:456)^(? + /i;i)^ 



(q + kiY (p + k456y 



Topology M: 



(6) , {p-ki2nq-k,,ycf^ + {p-k^)\q-kercf^ 

{p + k,,Y{q + k^,,Y 
, 4f ^ + {P- feO^cjf 4g + (g - A:e)^cr) , eg + (p - fcQ^c^^) 

The equations for topologies {B), (C) and (G) turn out to be redundant for the MHV 
amplitude with the choice of basis described above. 



B. Analytic Results for a 2 x 2 Block 

The structure of the equations is sufficiently complicated that it is not clear whether 
it is possible to find useful analytic solutions, so in practice we resort to solving them 
numerically. However the equations from topologies {D) and (F) are exceptionally simple 
and only involve the coefficients C2, c% and Cig, so they can easily be solved analytically as 
we now demonstrate. 
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1. Topology (F) 



The four contour integrals for topology (F) localize the integral at the four points 

(pW,g(-'))ij=i,2 given by 



Equation (13.51) is then 

(2) (2) 

„(6) , c^9 _ 4 „2 „ „(6) , ^19 _ n ('^^'^\ 

^2 + / fi) , , — '*'Si2'561 , C2 + t (2) , 1 ^2 ~ ^ ■ [6.16 

Note that even though there are four different contours we only obtain two independent 
equations since (13. 5p is independent of q. This is a generic feature whenever a contour is 



such that it chops off a massless box [16j. 
Solving (13.131) yields the coefficients 



(6) _ 4Si2S61 J2) _ 4Si2S61 . (1) , , .2 



C2 — , ^5 Cj^g 



(p(') + A:456)' (3.14) 



1 

where 

(p'" + fa«)^ (3.15) 



It is frequently useful to separate coefficients into their parity-even and parity-odd parts. 
Since parity exchanges (i, j) ^ j], it evidently takes a ^ 1/a. If we denote the parity 
conjugate of a coefficient c by c then we see that the even and odd parts of c'^' c^g'' are 
simply 

l(cf + cf ) = 2s,A , i(cf - cf) = 2».e»L (i±^ 

i(cS' + cg')=0, \[^S-^S)=i4.s yW'f . (3.16) 
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via 



The parity-even parts of these coefficients agree precisely with those obtained in 
the unitarity based method. Here we see that these coefficients can be obtained simply by 
solving two equations in two variables, and moreover the parity-odd parts automatically 
come along for free. 
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2. Topology (D) 



For topology (D) the contour integrals localize the integral at the points 

[2,6J 

= A6(aA6 + 7A3) q^'^ = T^AeAi , (3.17) 

(2,6) 



where 

_ £i3_+S23 ^ _ -S12 + a{si2 + S26) _ Su + ajsu + S26) fqisl 

""^36'^" (3|1 + 2|6] ' ^" [3|1 + 2|6) • ^^-'^^ 

Equation (IS.SI) then gives 

(6) (6) 

^19 + (p(l) _ fc^)2 = 4S12(S345S456 " S12S45) C^g^ + ^^(2) _ j^^y = • (3-19) 



We can eliminate cj^g to solve for Cg , finding 



^(4*^^ + ^) = 2si2(si23S345 " S12S45) 
1 /^(6) -(6)x ^ / X /C + 1 



-(Cg — Cg ) — 2Si2(Si23S345 — S12S45) ^— — -yj , (3.20) 

where 

_ (p(i) - kif 

Again the even part of Cg agrees precisely with the unitary based calculation of [l8|. Of 
course the equation fl3.19p provides a consistency condition on the coefficient c\q that we 
already obtained in (13.141) . 



C. The Parity-Even Part 



In the previous subsection we solved for the coefficients C2, Cg and Cig analytically and 



demonstrated that their parity-even parts agree with the results of 



18| . In order to check the 



validity of the leading singularity method it is important for us to compare the parity-even 



parts of all remaining coefficients as well. The first obstacle is the fact that [18| used a basis 
containing the integral shown in fig. HI According to our criteria we do not consider this a 
geometric integral since the propagator does not serve to cancel any pole. (The motivation 



for using / 



11 in H 



was that the integral is manifestly dual conformally invariant 



29| 
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hi (12) 



FIG. 4: An extra 12 integrals we add to the set of geometric integrals in order to facilitate com- 
parison with [isl]. These integrals are dual conformally invariant but not geometric. Nevertheless 
they can be expressed as linear combinations of the geometric integrals in fig. [3] using first class 
identities (see section IV). 

We show in the next section that reduction identities can be used to express In as a linear 
combination of geometric integrals, so secretly fig. [3] already contains In. However in order 
to facilicate comparison of our results with those of 18j it is convenient to explicitly add In 
to the set of integrals. Then we have a set of 177 + 12 = 189 integrals which is overcomplete 
by 18 + 12 = 30 elements. Encouragingly, we find that it is possible to find a solution of the 
equations in which the non-dual conformally invariant integrals /19-/24 all have coefficients 
whose parity-even part vanishes. This is a necessary condition for agreement with jQ] since 
the parity-even part of the amplitude was expressed there in terms of dual conformally 
invariant integrals. Moreover we find that after choosing this 'gauge' there is no further 
ambiguity in the basis; the linear equations furnish a unique solution. 

As indicated above, the equations are sufficiently complicated that we found it necessary 
to solve them numerically. Let us note however that by 'numerically' we always mean that 
we choose all of the spinors Aj, Aj to be random rational numbers (subject to momentum 
conservation, of course). Then the coefficients obtained by solving the equations always 
come out to be rational numbers, so they can be compared to the results of |3] with 
absolute precision. By repeated successful comparison for many different random values 
of the spinors, we are able to report complete agreement with the parity-even parts of 
the coefficients obtained from the leading singularity method with those obtained in 18| . 
recorded here in the (1) permutation: 

Cl = S61S34S123S345 + 5125455234-3345 + "5345(523-356 ~ ■S123S234) , 
C2 = 2512523 , 



C3 — -S234(-Sl23-S234 ~ -523-^56) 
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C4 


2 

— S12S234 , 




C5 


= -534(5123^234 " 


- 2S23S56) , 


C6 


= ~'Sl2'523'5234 5 




c? 


= 2S123S234S345 


- 4S61S34S123 - S12S45S234 - S23S56S345 


Cs 


= 2^61 (S234 -5345 


- S61'534) , 


Cg 


= "523'S34'S234 5 




ClO 


— i23l^''61''34 


''234''345j 5 


Cll 


= -5 12 -523 5234 , 




Cl2 


= -3345 (s 234 S 345 


— ■S61S34) , 


Cl3 


2 

— ■5345-^56 • 





However, as mentioned in the introduction we note that the leading singularity method is 
completely blind to the integrals /14 and J15 in [l^ since they have integrands that vanish 
in D = 4. 



IV. REDUCTIONS 



A. Methods of Reduction 



Repeatedly throughout the paper we have mentioned and used the new feature that hap- 
pens for six or more particles; loop integrals often satisfy linear relations which can be used to 



write one in terms of others (see for example 



20| , as well as 



for dimensionally regulated 



versions). Interestingly, there are relations even among what we call geometric integrals. 
Also important for us will be relations among integrals that appear in the manifestly dual 



conformal invariant expression of [18| and our geometric integrals. 

In this section we discuss in detail how these relations arise since it is a crucial step 
in completing the proof that the leading singularity method does determine the amplitude 
uniquely. Recall that out of the 177 coefficients the leading singularity fixes 159 thus leaving 
18 free parameters. We will now account for these as a consequence of relations or what 
we call reduction identities. In other words, the set of geometric integrals which naturally 
appears in the process of matching leading singularities is overcomplete. 

We distinguish between two different kind of identities; the ones that are valid on any 
contour of integration and the ones that are valid only on contours where loop momenta 
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can be taken to be in four-dimensions. We call these first and second class identities, 
respectively. 



1. First Class 

Identities of the first class are those which are valid on any contour of integration, in 
particular, on the real contours where integrals must be regulated, e.g. using dimensional 
regularization. 

First consider for example the pentagon-box integral 



with numerator factor {p — ki^, which is a permutation of the one shown in fig. HI Clearly, 
this integral is not geometric since the numerator is not a zero which cancels a pole of the 
integral. However, this integral is dual conformal invariant and it appears in the represen- 
tation of the even part of the amplitude obtained in |l8|. The goal is to write this integral 
as a linear combination of the geometric integrals in fig. [31 

Let us write the numerator as — 2ki ■ p. The first term cancels a propagator and gives 
rise to a double-box integral of type In in fig. [3l The second term can be decomposed 
by using that external momenta are kept in four dimensions. This means that only the 
four-dimensional component, p[4], of p contributes, i.e. ki ■ p = ki ■ p[4]. Given any four 
dimensional vector, one can write it as a linear combination of four axial vectors 



where t, with i = 1, ... ,4, form a basis of four-vectors. In the case of interest, we choose 
to write p[4\ as 




1/ 



\ 



^^lupa pi+1 pi+2 pi+3 



(4.1) 




(4.2) 



with the choice 



£i = ke, ^'2 = ^5 + ^6, h = ki + k^ + ke , £4 = -ki - k2 
This a standard construction in the scattering amplitude literature |20l. |28|. 
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(4.3) 



Since all of the £i are written in terms of external particle momenta, they are completely 
four-dimensional so we are free to replace p[4] by p in the coefficients (p[4] • ii). Writing 
each coefficient as 2p ■ £i = {p + ii)'^ — p^ — £f , one finds that the term fci • p in the original 
integral can be decomposed in terms of numerators which give rise to geometric integrals. 
Perhaps the only term which might require some explanation is the one corresponding to 
£4 = —ki — k2 since it is one which does not cancel a propagator. In this case the factor 
(p — ki — k2f' becomes a numerator which is easily seen to be a zero that cancels a pole 
which removes a leading singularity and hence gives rise to a geometric integral. 

Prom this example it is clear what the necessary conditions are for the existence of first 
class relations among integrals. The first condition is that there be at least four propagators 
(including involving the same loop variable and only external momenta. If the number 

of such propagators is exactly four then there must be at least two different ways of putting 
a numerator which only involves the loop variable of interest and external momenta. This 
was the case considered above. If the number of propagators is at least five, then a relation 
can be obtained if at least one numerator is available. If the number of propagators is six 
or more, then no numerators are needed. 

In our case, with six external particles, the maximum number of propagators in a single 
loop which only depend on a single loop variable and external momenta is five. This diagram 
is a hexagon-box. In fact, it is easy to show that a hexagon-box with a numerator can be 
reduced. We leave this as an exercise for the reader since such an integral did not have to 
be included in the original set of geometric integrals, at least in the MHV case. 

Using identities of the first class we will be able to show that all six /12 integrals can 
be written in terms of /13 integrals plus other geometric integrals. This shows 6 of the 18 
relations we have to account for. Also using first class relations we will show that 6 of the 
12 integrals can be written in terms of the remaining six and other geometric integrals. 
Summarizing, after using all first class identities we end up with 6 relations left to explain. 
These turn out to be second class identities. 

2. Second Class 

The identities discussed above are valid in any contour because the loop momentum is 
always contracted with external momenta and hence it can be treated effectively as four 
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dimensional. We now turn to identities which only hold on the contours where the 
loop momentum integrals can be taken completely in four dimensions. These identities will 
not hold in dimensional regularization. In fact, their failure to hold exactly is precisely 
proportional to integrals with numerators made out of the — 2e-dimensional component of 
the loop momenta. We obtain identities by reducing integrals of this type to find linear 
combinations of other integrals which must sum to zero for four dimensional loop momenta. 

Consider first an integral with numerator, P[-2e] ■ ^[-2e] = V ' 1 ~ P[4] ' ■ The first term 
can be written as, 2p ■ q = {p + qY — — . In order to treat the second term one has to 
find a convenient basis to expand p[4] and one for qy^^ . The relevant diagram must have at 
least four propagators that only depend on p and at least four that only depend on q and 
external momenta. With six external particles there is only one possible diagram (up to 
relabeling). This is the double pentagon integral 




Let us choose to write pj^j using the reference vectors 

£i = fee , ^2 = ^5 + ^6, l^ = ki + h + h, £4 = ki (4.4) 
while writing gj^j using 

£1 = fci , 4 = fci + ^2 , i3 = ki + k2 + ks, £4 = kQ. (4.5) 

Plugging these into (14. 2p and calculating j9[4] ■ g[4], we find an expansion containing ge- 
ometric integrals and some further integrals that can very easily be further decomposed 
into geometric ones. These identities give rise to six relations among the remaining six I24, 
integrals and the rest of the geometric integrals. The reader might wonder how can one get 
six equations if the starting point, which is the P[-2e] ■ 1l-2e] pentagon-pentagon integral, has 
a 4-fold symmetry implying that there are only three independent such integrals. However, 
the decomposition process breaks that symmetry since one needs to make a choice of basis 
vectors as in equations (14.41) and (14. 5p . There are two independent choices, leading to a total 
of six independent reduction identities. 

One could also consider integrals with a factor of pf_2e] = — P^4] in the numerator. 
However the only way to use a reduction of such an integral in the case of six particles is 
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if p is the loop momentum in a hexagon inside the hexagon-box integral. The hexagon-box 
does not appear in fig. [3] since it is never needed in order to solve the equations, so we have 
no need for such reduction identities. 

B. Summary of All Integral Reductions 

Here we summarize all of the relevant reduction identities that can be derived using the 
techniques explained above. First we have the identity schematically represented as 




which we use to indicate that the integral on the left can be expressed as a linear combination 
of the integrals on the right. It is straightforward to work out all of the precise coefficients, 
but they are not important for our analysis. The important point is the conclusion that the 
integral shown in fig. H] can be reduced to integrals already present in fig. [31 



It is also straightforward to derive the first class relation 




This is again a schematic relation: the integrals on the right-hand side can appear in various 
rotated or fiipped incarnations. 
The final first class relation is 




which implies that of the 12 apparently independent integrals of the type J24 shown in fig. [3l 
in fact only 6 are linearly independent. 

Next we summarize the second class reduction formula, discussed in section IV.A.2, that 
is obtained by starting with a double pentagon with P[-2t] ■ ?[-2e] numerator. As explained 
above this analysis leads to 6 independent identities. Schematically these identities take the 
form 




Again the integrals on the right-hand side can appear in various different permutations. 

Talcing everything into account, we find that the set of geometric integrals in fig. [3] is 
overcomplete by 6 + 6 + 6 = 18 elements. This precisely accounts for the 18 free parameters 
we found in solving the linear equations for the MHV coefficients. (This increases from 18 
to 18 + 12 = 30, as discussed in section III.C, when the integral In is thrown into the mix.) 



V. NUMERICAL CHECK OF THE ABDK/BDS ANSATZ 



One of the most interesting properties of scattering amplitudes in Ai^ = 4 SYM is that the 
structure of infrared divergences in higher loop amplitudes is very simply related to those 



of lower loop amplitudes 



32l |. In 22|, |23|, it was conjectured that this simplicity persists, at 



least for MHV amplitudes, to the finite terms as well. The precise form of the conjecture at 
two loops, in dimensional regularization to D = A — 2e, is 



(2) 



(1) 



n,MHV 



(e)2 + /(2)(,)M, 



4iHv(2e) - ^ + 0(e) 



where M 



(5.1) 
-(C(2) + 

C(3)e + C(4)e^ + ■ ■ ■)■ The conventions implicit in equation (15. ip require that every loop 
momentum integral p be normalized with the factor 



n.MHV 



n.MHVV^ 

^i^MHvMn.MHV the normalized L-loop amplitude and /'■^-'(e) 



MHV 



(5.2) 



The simple structure (15.11) holds perfectly for n = A and n = 5 1 171 . 1221 . |2J]. However 
it apparently fails beginning at n = 6 particles. This was found in jlSj] by computing the 

(2) 

parity-even part of Mg ^-^y numerically and finding disagreement with the right-hand side 
of ( 15. ip . Here we do not have anything to add to this issue except that the parity-even piece 
of our full answer agrees with the result of 3], thus providing independent confirmation. 

Since the leading singularity method allows us to obtain the parity-odd parts of all co- 
efficients with no more effort than the parity-even parts, we are in a position to test the 
parity-odd part of (15. ip for n = 6. Note that the loop momentum integrals must necessarily 
be evaluated numerically (except for n = 4, where analytic results are known through three 
loops) with current state-of-the-art technology (in particular we use 

Restricting to the parity-odd part of ( 15. ip yields 
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3l|). 



(2) 



6,MHV odd 



(e) = M, 



(1) 

6,MHV even 



;e)M, 



6,MHV odd' 



e) + 0{e). 



(5.3) 
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The one-loop amplitude is [Gj 



r(l) 

'6,MHV evenV^; — 



i=l 

6 



-^6,MHV odd(^) 



~\ Yl ((^K+lK+2K+3|i+4] - [z|z+l|i+2|z+3|i+4))Pi+5,i+6(e) ,(5.4) 

where Pj+5,j+6(e) is the one- loop pentagon integral with external legs i+5 and i+6 joined 
and with a factor of pf_2e] ^^e numerator. 

We have evaluated (15. 3p numerically at two independent kinematic points with randomly 
generated values of Aj and Aj for the six external particles. Denoting the left- and right-hand 
sides of (15.31) by L and R respectively, we find 

^ _4xj^ ^ 4x^ ^ 1(2) X 10- _ 0430(7) _ ^ 

Rie) = _ ^^i^ _ 0.92(1) + 0(e) (5.5) 

at the first point and 

i(,) = 1^ _ X f " + _ 15.902(5) _ ^ ^ ^ _ 

g4 g2 g 

i?(e) = - IMM^ _ 60.38(4) + 0{e) (5.6) 

at the second. In these expressions the value in parentheses denotes the estimated numerical 



error in the last digit as reported by 



3l|. 



We emphasize that the cancellation of the divergent terms in (15. 5p is not a check of the 
ABDK/BDS conjecture, but rather a check on our application of the leading singularity 



method. This is so because it is a known fact 32], not a conjecture, that (15. ip must hold 
for the infrared divergent terms. Had we gotten a nonzero result, it would have signalled 
an error in our calculation of the integral coefficients. Note that the cancellation in L(e) 
is highly nontrivial in the sense that the result shown is obtained after summing of order 
100 contributions which are typically of order 1. It is the fact that we see the cancellation 
persisting to order 0{e^) that strongly suggests that the parity-odd part of the two-loop 
six-particle MHV amplitude indeed satisfies the ABDK/BDS conjecture (15. 3p . 

However it is important to note once again that since the leading singularity method is not 
sensitive to any "/i" -terms we cannot rule out the possibility that there may be additional 
such contributions to Mg^^inv odd- Given our apparently successful check of (15. 3p there are 
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three possibilities: (1) the two-loop amplitude does not contain any parity-odd /i-terms (this 
is indeed the case for = 5 particles [17]), (2) the amplitude does contain parity-odd /x- 
terms but they contribute only at C(e), or (3) the amplitude contains parity-odd /i-terms 
which spoil the ABDK/BDS relation flETD . 

VI. CONCLUSION 

In this paper we employed the leading singularity method to determine the integral 
coefficients of the planar two-loop six-particle MHV amplitude in A/" = 4 YM, the parity- 
even parts of which were recently obtained in 1^ using the unitarity based method. One 
advantage of the leading singularity method is that the full coefficients, including parity-odd 
parts, emerge from solving the relatively simple set of linear equations displayed explicitly 
in section III. A. 



The leading singularity method has previously proven succesful 16| at reproducing the 



n = 4 33| and n = 5 17|] particle amplitudes at two loops. However, there is currently no 
proof that a general amplitude in A/" = 4 Yang-Mills is uniquely determined by its leading 
singularities only. It is a logical possibility that finding a representation of an amplitude in 
terms of simpler integrals which faithfully reproduce all of the leading singularities is not a 
sufficient condition to guarantee correctness of the representation, although clearly it is a 
necessary condition. 

Here we find that the n = 6 particle amplitude is in fact completely determined by its 
leading singularities, although establishing this fact required that considerable attention be 
given to the choice of basis for the integrals and reduction identities which relate various 
integrals to each other. This was necessitated by the fact that the full set of linear equations 
we found does not have a unique solution. Fortunately we found that all of the ambiguities 
could be accounted for by taking into account reduction identities. 

One nice feature of the leading singularity method is that the procedure naturally provides 
a set of 'geometric' integrals for any amplitude under consideration. The set of geometric 



integrals does not coincide with the manifestly dual conformally invariant [27|, |29|] basis used 
to express the parity-even part of the n = 6 amplitude in jl8|. We expect this to be true in 
general. This is not surprising given the fact that by using the leading singularity technique 
both even and odd parts of the amplitude are computed simultaneously, whereas already for 
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five particles the odd part of tlie amplitude 17|] is not expressible in terms of integrals with 
manifest dual conformal properties. Our results indicate that also for n = Q the odd part of 
the amplitude cannot be expressed in terms of dual conformally invariant integrals alone. 

Another motivation for computing the MHV six-particle amplitude, beyond its serving 
as a testing ground for the leading singularity method, is to study the so-called ABDK/BDS 



conjecture for MHV amplitudes 
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231] . Although the conjecture was shown to be violated 



by the parity-even par t of the n = 6 particle amplitude jlSj], following earlier doubts that 
had been raised in 



34 



35 



36j, we provide numerical evidence that the parity-odd part of 



the amplitude does satisfy the ABDK/BDS relation. 

Equivalently, one can say that the parity-odd part evidently cancels out when one takes 
the logarithm of the resummed amplitude (see 23|). Although we do not know of any proof 
that this has to be the case, the result is consistent with the structure seen at strong cou- 



pling 



37l |. which is manifestly parity- invariant. It is also consistent with the astounding 



but still mysterious equivalence between scattering amplitudes and lightlike Wilson loops in 
A/" = 4 YM that has been observed at one-loop 29|, [381] and at two-loops through n = 6 par- 



ticles 39, 
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42], since the 'vanilla' Wilson loop does not carry any helicity information. 
It is a very interesting open problem to determine if the amplitude/ Wilson loop equivalence 
can be extended to other helicity configurations by appropriately dressing the Wilson loop. 
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